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Y ■ Abstract 

fH , The phase ordering dynamics of coupled chaotic bistable maps on lattices with defects 

Q ' is investigated. The statistical properties of the system are characterized by means of 

O ' the average normalized size of spatial domains of equivalent spin variables that define the 

phases. It is found that spatial defects can induce the formation of domains in bistable 
I ' spatiotemporal systems. The minimum distance between defects acts as parameter for a 

^ \ transition from a homogeneous state to a heterogeneous regime where two phases coexist 

■ The critical exponent of this transition also exhibits a transition when the coupling is 
' increased, indicating the presence of a new class of domain where both phases coexist 
, forming a chessboard pattern. 

(N 

' Coupled map lattices constitute fruitful and computationally efficient models for the study 
. of a variety of dynamical processes in spatially distributed systems [1]. The discrete-space 

■ character of coupled map systems makes them specially appropriate for the investigation of 
I spatiotemporal dynamics on nonuniform and on complex networks [2-5]. 

' There has been recent interest in the study of the phase-ordering properties of systems 
I . of coupled chaotic maps and their relationship with Ising models in statistical physics [6- 
'"O ' 12]. These works have generally assumed the phase competition dynamics taking place on a 
^ . uniform space; however, in many physical situations the medium that supports the dynamics 
can be nonuniform due to the intrinsic heterogeneous nature of the substratum such as porous 
or fractured media, or it may arise from random fluctuations in the medium. This paper 
investigates the process of phase ordering in coupled chaotic maps on a lattice with defects as 
a model for studying this phenomenon on nonuniform media. 
^ ■ We consider a system of coupled maps defined on a two-dimensional square lattice of size 
N ~ L X L with periodic boundary conditions and having randomly distributed defects, as 
shown in Fig. 1. A defect is a non-active site, i.e., a site that possesses no dynamics. The 
density of defects in the lattice is characterized in terms of the minimum Euclidean distance d 
between defects. The dynamics of the diffusively coupled map system is described by 

x.(t + l) = (l-e)/(a;,(t))-f -^^/(x. W) , (0.1) 

where Xi{t) is the state of an active site i (i = 1, . . . , N) at time t, ui is the set of the nearest 
neighbors of site i and Mi G {1,2,3,4} is the cardinality of this set, e measures the coupling 
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strength, and f{x{t)) is a chaotic map that expresses the local bistable dynamics [7,8], 

(-2^/3 -fix, if xe [-1,-1/3] 
fix)=l fix, if X £[-1/3,1/3] (0.2) 

[ 2^/3 -Aia:, if xe[l/3, 1]. 

For fi G (1,2) the map has two symmetric chaotic band attractors, one with values Xi{t) > 
and the other with Xi{t) < 0, separated by a gap about the origin. Then the local states have 
two well defined symmetric phases that can be characterized by spin variables defined as the 
sign of the state at time t, ai{t) — sign{xi{t)). 




Fig. 1. (Left) Spatial support of the system, showing active sites (o) and defects (•). Defects are 
randomly placed in such a way that their density distribution is maximum for a given minimum 
distance d between defects. (Right) The defect density scales as p = 0.625d~^. 

We fix the local parameter at /i = 1.9 and set the initial conditions as follows: one half of 
the active sites are randomly chosen and assigned random values uniformly distributed on the 
positive attractor while the other half are similarly assigned values on the negative attractor. 
If the number of active sites is odd, then the state of the remaining site is assigned at random 
on either attractor. 

In regular lattices {p = 0) phase growth occurs for values e > eg, where — 0.67 [8]. In 
contrast, in the defective medium there exist a minimum value of p for which the domains 
formed by the two phases reach a frozen configuration for all values of the coupling e. To 
characterize the phase ordering properties of the system Eq. (10. 1|) we use the normalized size 
of the phase domains, averaged over 50 realizations, as an order parameter given by 

r = l ij 

where rij is the Euclidean distance between nodes i and j. Figure [5] (left) shows the variation 
of R in the space of parameters (e,d). Note that the system is heterogeneous (i? ^ 0), with 
coexistence of the two phases for values e < eo. When e > cq two regions appear in the 
phase diagram: an heterogeneous regime for small values of d (large p) and an single-phase, 
homogeneous state {R = 1), for large values of d (small p). 

The transition between these two types of behaviors induced by the presence of defects for 
a given value of e > eo can be described by the scahng relation i? ~ (dc — c?)^^'\ where dc{e) is 
the threshold value of the minimum distance at which the transition occurs, and /3 is a critical 
exponent. 

Figure 2 (right) shows how the behavior of the critical exponent /? varies with e. There 
is a critical value of the coupling Cc ~ 0.855 at which (3 changes from being a constant value 
of /3 w 0.24 to follow a scaling relation /? ~ (e — ec)'', where 7 « 0.25 is a critical exponent 
characterizing the transition. This transition is related to the emergence of a new type of 
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Fig. 2. (Left) Normalized domain size Ji as a function of e and d. Dark color represents homogeneous 

(single-phase) states and light color corresponds to heterogeneous (two-phase) states. (Right) Critical 
exponent (5 vs. coupling e. The two inserts show typical frozen configurations of the system for e < tc 
(left) and for e> tc (right). Dark and light colors represent each phase and black dots correspond to 
defects. 



configuration of domains induced by defects in the system, where coexistence of alternating 
sites of the two phases takes place forming a chessboard pattern. It should be noticed that 
continuous variation of the values of critical exponents may also occur in some statistical 
mechanics models showing phase separation dynamics [13, 14]. 

In summary, we have found that spatial defects can induce the formation of domains in 
bistable spatiotemporal systems. The minimum distance between defects acts as a parameter 
for the transition from a homogeneous to a heterogeneous regime. The critical exponent of 
this transition (3 exhibits a second order transition when the coupling is increased, indicating 
the presence of a new class of domain where both phases alternatively coexist. 
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UNET, F-2002000426 from FONACIT, and 1-886-05-02- A from CDCHT of Universidad de Los 
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